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Abstract 



In paper I, we have studied the conditions under which flux tubes are expelled 
from adiabatic thick accretion disks. In the present paper, we explore a few other 
models of thick disks, where flux tubes could be stored. We show that flux tubes 
with sufficiently weak fields are not expelled out if they move adiabatically inside 
an isothermal disk; they continue to oscillate around mean equipotential surfaces 
inside the disk. If the field in the flux tube is amplified due to the shear, they 
are eventually expelled away. We explore a 'toy' model also, where the entropy 
increase outwards from the center of the thick disk and find a similar behavior. 
Flux storage in the disk, as in the case of the sun, in general, enhances the 
possibility of sustained magnetic activity formation of coronae in the chimney 
region. The existence of coronae on the disk surface may explain the short- 
time variability in the spectra of Blazars and the emission of energetic particles 
from AGNs and Quasars. It may also supply matter to the cosmic jets through 
magnetized winds. 
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1. INTRODUCTION 

In an earlier paper (Chakrabarti & D'Silva 1993; hereafter referred to as Paper 
I), we have presented a detailed study of how magnetic flux tubes, either generated 
inside a disk or advected along the accreting flow, are expelled outside the disk. Over 
a large region of the parameter space, spanned by the size and the fleld strength, the 
flux tubes are seen to emerge in the chimney region of the thick accretion disk. These 
flux tubes emerge in buoyancy timescales, and might be advected along the jets, which 
are beheved to be originated in the chimney of the disk. It is long known that most of 
the luminosity of the disk is emitted from the chimney. For the first time, wc pointed 
out that the chimneys could also be magnetically the most active regions in the disk. 

Blazars are long believed to belong to that class of active galaxies where the jets 
point at the observers. These objects are known to be highly variable over the entire 
electromagnetic spectrum. Variability time scales vary from lO^s to lO^s (e.g., Bregman 
1988). A large number of models have been proposed to explain the cause of variability 
(see. Paper I for a review) . If it is caused by the magnetic activity at the base of radio 
jets (Blandford and Konigl 1979), then it is imaginable that the flux tubes are actually 
anchored in the chimney region of the thick disks and the activities are due to shock 
production by flaring events. 

In the case of the sun, flux tubes are believed to be anchored beneath the convec- 
tion zone where the temperature gradient is sub-adiabatic, and sunspots are believed 
to be formed at the regions where these flux tubes emerge from the surface (Parker 
1979). This emergence takes place due to Parker's instability (Parker 1955), which is 
the common cause of the formation of sunspots, coronae, prominences, flares and other 
magnetic regions on the solar surface. In this paper, we show that thick accretion disks 
could also entertain a very similar scenario where flux tubes can be anchored inside 
the disk and the parts of these flux tubes could emerge into the chimney, making it 
a region of sustained magnetic activity. So far no work has been done to show that 
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flux tubes could be anchored inside a thick accretion disk. In the present paper, we 
explore models of thick disks which render efficient anchoring. We find that isothermal 
disks are capable of storing flux tubes provided the tubes move adiabatically. We also 
explore some 'toy' models of thick disks, which are in hydrostatic equilibrium, and at 
the same time, have entropy rising outward from the center. These disks can store flux 
tubes provided the field strength is sufficiently weak so that even the buoyancy force 
is counter acted. In other words, in order to anchor a flux tube of given strength the 
local entropy gradient should be sufficiently sub-adiabatic. In may be noted in this 
context that Shibata et al. (1989), using a two dimensional MHD simulation show 
that the Parker instability of horizontal magnetic flux sheet could be suppressed inside 
a two temperature inversion medium. Such models could be physically realised in a 
variety of astrophysical circumstances which we elaborate. Axisymmetric flux tubes, 
in the absence of shear ampliflcation and other dynamo effects, continue to oscillate 
around mean equipotential surfaces indefinitely, without ever being expelled out of the 
disk. Perturbative effects, such as deviation from axisymmetry (either of the fields, or 
temperature or pressure profile) can easily bring these tubes to the surface and cause 
magnetic activity as in the case of the sun (Moreno-Insertis, 1986; Choudhuri, 1990; 
D'Silva & Choudhuri, 1993). In the context of accretion disks, such perturbations 
could be caused by various interchange instabilities (Kaisig, Tajima & Lovelace, 1992) 
or non-axisymmetric shear instabilities (Balbus & Hawley, 1991, 1992; Hawley & Bal- 
bus, 1992). Shear amplification of the field could also enhance the field strength till it is 
expelled away due to Parker instability (e.g., Matsumoto et al. 1990; Shibata, Tajima, 
Matsumoto 1990). 

The 'toy' models we explore have a simple entropy distribution, wherein the entire 
disk becomes radiative and fiux tubes can be stored, or the entire disk is convective 
and all fiux tubes get expelled out of the disk. In principle, one could try out entropy 
distributions where the disk can have a radiative core and a convective envelope, as 
in the sun. This exercise leads to the most obvious question - 'What governs the en- 
tropy distribution in an accretion disk?'. The most serious problem is our ignorance 
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of viscosity and its role in controlling the dynamics of the entire disk. For instance, 
if accretion was surface phenomenon (Paczinski 1978) in thick disks, then, naively one 
would expect that the ratio of gas pressure to radiation pressure (5 would be small at 
the surface and increase inward to the center. Since the entropy function K (described 
in Paper I) of thick disks goes roughly as for /3 << 1, this allows the entropy to 

rise from the center to the outer edge of the disk. As we shall show later, this is the 
Schwarzchild criterion for convective stability in thick disks. On the other hand, if ac- 
cretion is equatorial (Paczyhski and Wiita 1980; Wiita 1982, Abramowicz & Paczyhski 
1982, Chakrabarti, Jin & Arnett 1987), the disk can naively be expected to be convec- 
tive. Convective stability also implies stability towards magnetic buoyancy, provided 
the field strength is low enough — the value depends on the entropy gradient. 

The presence of magnetically active coronae on disk surface can not only explain 
the variability of Blazars, they can also supply matter to the jets through magnetic 
winds. The flares can accelerate particles through the first order Fermi process as in 
the case of the sun. Thus a large fraction of the high energy radiation from AGNs may 
originate at these flares. Magnetic flux tubes also preferentially fllter charged particles 
(e.g. protons, electrons and ions) from those uncharged (e.g., neutron) and may produce 
neutron rich core of the tori which, in turn, can produce neutron rich elements through 
processes involving rapid capture of neutrons. 

The plan of the present paper is the following: In the next Section, we describe 
model equations for the thick disks which we study in this paper. The equations of 
the dynamics of flux tubes inside a disk have already been presented in Paper 1. Here 
we just reproduce the equations and show how the quantities inside the flux tubes 
vary, when the tubes move isothermally or adiabatically. In §3, we persent numerical 
solutions describing the behavior of flux tubes in the various thick disks. In §4, we briefly 
describe the astrophysical implications of our results. Finally, in §5, we summarize our 
paper and make concluding remarks. 
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2. MATHEMATICAL EQUATIONS 

We consider axisymmetric flux rings, and study their dynamics under the thin 
flux tube approximation, both of which are described in Paper I. Here we reproduce 
Eqns. (10)-(12) of Paper I. 



and 
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^9 + 2^9 + — [-^02 sin 9 cos 9- 2^u^ sin 9 cos 9] = 
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= - 7T^Vt{ ^ cos ^ + Tens cos 9 + } , (2) 
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where (^, ^, 0) describes the position of a fluid element inside the flux tube, and the 
quantities of drag D, dimensionless angular velocity oo and X are described in § 2.2 of 
Paper I. Magnetic tension is given by 



T = ^-^^ (4) 



The expansion factor A — (a/aoY is the ratio of the cross-sectional area of the flux 
tube with its initial cross- sectional area, and can be derived using the mass and flux 
conservations as 
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. ^ Pe(^o,^o) . ^o singp (I-Mq) 

The subscripts e and i refer to the corresponding quantities outside (ambient medium) 
and inside the flux tube respectively. The subscript o refers to the quantities corre- 
sponding to the initial conditions. 

2.1 Isothermal thick disks 

In the case when efficient cooling mechanism for ions is not possible, ions inside 
the disk may remain very hot, close to the virial temperature. Such disks have been 
invoked in the context of two temperature tori (Lightman and Eardley, 1974) and ion 
pressure supported tori (Rees et al. 1982). For the purpose of the present analysis, 
we consider isothermal disks, as another extreme from the adiabatic case discussed in 
Paper I. Gas pressure obeys the equation of state, Pg = RpTd/ fi, Td being the constant 
temperature of the disk. The radiation pressure Pr — {l/3)aT^ remains a constant. 
The equipotential surfaces follow the same equation as for adiabatic disks (Eqn. 2 of 
Paper I). The density variation of the disk obeys the distribution, 

pK) ^ p(fljexp{-^i:^(-l^ - ^) + _ (J.i„,)--1), (6) 

The flux tubes may radiate very efficiently and remain in isothermal equilibrium 
with its local surroundings. On the other extreme, when the cooling is not efficient, 
the flux tubes can move adiabatically inside the disk without exchanging any energy. 
Sections 2.1.1 and 2.1.2 give the relevant equations for these two cases. Throughout the 
calculations we choose an initial entropic condition in which the flux tube is in thermal 
equilibrium with its surroundings. The condition of pressure equilibrium provides the 
initial buoyancy factor to be Mq = B'^/ {Sirpe) ■ 

2.1.1 When the flux tubes move isothermally 

We now consider the case where the flux tubes are in thermal equilibrium with 
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its surroundings. This is possible when the exchange of energy is very efficient. The 
radiation pressures inside and outside the tube are equal. Using the mass and the flux 
conservation {B/pi^ sin 9 — constant) and the condition of pressure equilibrium, we get 
a quadratic equation in M, whose solution is of the form 



This equation is similar to Eqn. (15) of Paper I. The solution with a positive sign is 
unphysical. 

2.1.2 When the flux tubes move adiabatically 

In this case, the entropy of the flux tube remains a constant throughout its motion 
inside the disk. Using the mass and flux conservations along with the condition for 
pressure equilibrium, we get a polynomial in Pi/pe, 




(7) 



where 
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(9) 
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Pi 



;i-MoKo + (aAi/3i?)r3„ 
2.2 A toy model for the radiative thick accretion disk 

Here we present a model of a thick accretion disk in which the entropy is a 
function of the potential W of the disk. In the adiabatic disk, entropy is independent of 
W and in the isothermal disk S (x —log{p). Presently, we use a simple functional form, 
wherein, either the entire disk remains radiative or the entire disk becomes convective, 
depending on whether the Schwarzchild criterion for convective stability (Landau & 
Lifshitz 1959) is satisfied or not. In thick accretion disks the Schwarzchild criterion for 
convective stability can be written in a modified form as 



dK , , 



We take a linear form for the entropy function K as 



K[W) =K^ + K2W , (12) 



where Ki and K2 are constants. 
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(14) 



Wc = W{Rc, 9 = 90°) is the value of W at the center of the disk and Wont = W{Rout, = 
90°) is its value at the outer edge of the disk which is chosen to be Rout = 20rg. This 
dependence of K on W could be through jS or through fi or both. We assume that the 
composition of the disk is homogeneous. The local value of the entropy function K{W) 
is given by 



where l3{W) is the ratio of the gas pressure Pg to the total pressure pt on a given 
equipotential surface, Rg is the gas constant, a is the radiation density constant and 
fi is the mean molecular weight of the gas. According to the Schwarzschild criterion, 
if K2 > 0, the disk is radiative, and if K2 < 0, it is convective. K2 = gives the 
adiabatic disk. The potential W increases outwards from the center of the disk, so 
{Wc - Wout) < 0. In Eqn. (|l|), K2 > for K{Wout) > K{Wc); the disk is convectively 
stable (radiative disk). For a radiation pressure dominated disk [jS << 1) entropy goes 
roughly as /3~'^^^. By choosing (3{Wc) > (3{Wout) the disk becomes radiative {dK/dW > 
0), or by choosing P(Wc) < PiWout) it can be made convective {dK/dW < 0). 

As in Paper I, we choose the polytropic equation of state P = Kp^ inside the 
disk, with K = K{W). We assume that the disk is a radiation pressure dominated 
geometrically thick disk. It is also optically thick, photon and matter act as a single fluid 
of polytropic index 7 = 4/3. We also assume that the equation of state is barotropic, 
P = P{p) and the force balance equation (Eqn. 1 of Paper I) is integrable W{P) = 



K{W) = [ 



(Ijl^W) 3 M^ll/3 



(15) 
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— / (dP/p), and 



»'-»'o., = -^-^K.ta«)--, (16) 



and p is given by 



2.2.1. When the flux tubes move isothermally 

As usual we assume that the initial entropic condition is such that the flux tube 
is in thermal equilibrium with its surroundings. Thus, 



^ = 1 - Mo. 

If the flux tubes remain in thermal equilibrium throughout their journey in the disk, 
then the pressure balance condition with the mass and flux conservations give the 
expression of the buoyancy factor M to be the solution of a quadratic as in Eqn. (|^), 
where the constant h is 



lj = 2+'^° sin^„^2,Pe(eo,^^o),2(l-^o)'/5e,o(l-/3e 



esin^ ^ ^ p,{i,e) ' Mo (3eil-Pe,o) 

2.2.2. When the flux tubes move adiabatically 

In this case, Pi{^,0) = [3i{^o,6o). Using the mass and flux conservations we get 
an equation for Pi/ Pe in the similar form as Eqn. (|), where the constants ki and /c2 
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are given by 



1 - f3eAW)Mo Ke, 
(1 - 7V/o)4/3 



(18) 
(19) 



2.3 Inclusion of effects of shecir amplification 



The field strength of the flux tubes can be amplified by convective motions and dif- 
ferential rotation, and the azimuthal field can be shown to grow exponentially (Galeev, 
Rosner & Vaiana 1979), with a growth time r ~ {C,/v), where v is the convective ve- 
locity. Following Bisnovatyi-Kogan and Blinnikov (1977), we approximate convective 
velocities to be ^; a^^^Cg, where Cg is the sound speed and is roughly equal to the local 
infall velocity v ~ {ck/\^^ — 1); a being the viscosity parameter (Shakura & Sunyaev 
1973). 

Assuming that the shear amplification of field inside a thick disk takes place at 
the same rate as inside a thin disk, we choose. 



B = B'exp{t/T}, (20) 

where B' is the field in the previous time step and t is the time interval. It is not 
our intention to develop a detailed picture of the dynamo process inside the disk. Our 
present strategy is to get some insight about the physical processes that might take 
place when such an enhancement of the field is turned on. 

When the field is amplified, flux inside the tube is no longer conserved, so the 
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Eqns. @ and (|T^) change accordingly as 



ens 
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(21) 



and 



[l-/3e,o(iy)Mo]< 
(1 - Mo)4/3 



(22) 
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e \2 



(23) 



where, a prime (') refers to the corresponding quantities at the previous integration 



We first study the isothermal disks in § 3.1 and then the radiative and convective 
disks in § 3.2. 

3.1 Isothermal disks 

Since we are interested in studying the generic behavior of the flux tubes in 
isothermal disks, we choose some typical parameters. We consider an isothermal thin 
disk described by Eqn. (^), with the boundary condition of p{Rout) = 10^^^ gm cm~^. 
The disk is around a black hole of mass of 10^ M©, with its outer edge Rout at 20rg. 
The temperature of the disk is taken to he = 5 x 10^ °K. The angular momentum 
distribution / = /o(^ sin^)", where lo and n are positive constants as in Paper I. In § 3.1, 
we study the behavior of the flux tubes when they are in isothermal equilibrium with 



step. 



3. NUMERICAL RESULTS 
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their surroundings and § 3.2 gives the behavior when they move adiabatically inside 
the disk. 

3.1.1 When the flux tubes move isothermally 

We take a constant angular momentum disk {n = 0) with lo = 2. Flux rings of 
Mo = 0.3 are released close to the equatorial plane (6* = 89°) at ^ = 4, 5, 6, 7 and 10 with 
the initial conditions of zero velocity, and the trajectories of the flux tubes are computed 
by integrating Eqns. ([l|)-(0) with the assumption that drag and accretion are negligibly 
small ~ and Uacc ~ 0). Figure la shows the results; it is almost identical to the 
corresponding adiabatic case in Paper I, except that the flux tube released at ^ = 6 also 
emerges into the chimney. This interesting behavior is because magnetic tension (which 
goes inversely to the radius of the flux ring) dominates over magnetic buoyancy. Section 
3.3 of Paper I discusses the importance of magnetic tension. Simple analytic work which 
involves balancing magnetic tension and the buoyancy effects, shows that if > 4x lO^'^ 
°K, magnetic tension dominates over magnetic buoyancy and the flux tubes released at 
the outer edge of the disk should also fall into the chimney. Numerically, when all the 
effects are taken into account, this happens at somewhat lower temperature, even at 
Td = 10^° °K itself, as shown in Fig. lb. All the conditions in Fig. lb, except the value 
of Td, are the same as in Fig. la. 

In order that the self-gravity of the disk does not become important, we chose 
disks of very small density, i.e., we choose a very high temperature. The black hole 
mass Mbh = 10^ Mq and the disk of size of 20rg with the outer boundary condition 
of p{Ro) = 10^^^ gm cm~^ restrict the lower limit on the temperature of the disk to 
roughly 10^ °K. On the other hand, if is increased beyond 10^^ °K the disk becomes 
ion pressure dominated. The disk parameters of Mbh, p(-Ro) and -Ro that we have 
chosen, conspire in such a way that Td has to lie between 10^ and 10^° °K in order 
to avoid magnetic tension from dominating the dynamics and the disk becoming self- 
gravitating. A smaller temperature can be chosen provided we choose a smaller disk 
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size. When is below 10^° °K, the behavior of the flux tubes is similar to that of 
the tubes which move isothermally in the adiabatic disk (§ 3.1 of Paper I). However, 
when flux tubes move adiabatically in the isothermal disk the behaviour is remarkably 
different as discussed in the next subsection. 

3.1.2 When the flux tubes move adiabatically 

We repeat the calculations with parameters as in Fig. la, except that the flux 
tubes move adiabatically inside the disk. This is done by calculating the buoyancy 
factor M using Eqn. (^. To start with, the flux tubes are buoyant because of the 
particular initial condition chosen such that they are in thermal equilibrium with their 
local surroundings, but not in hydrostatic equilibrium. The flux tubes move out in the 
direction of the pressure gradient force and in doing so they expand in order to be in 
pressure equilibrium with its immediate surroundings. The reduction in density due to 
the expansion, however, is not sufficient to keep its density smaller than its surroundings, 
and the tube becomes less buoyant as it rises, till it reaches a point where the tube 
is no longer buoyant (pe = Pi)- The rising flux tube overshoots this point and enters 
into a region where it becomes heavier than its surroundings and begins to fall back 
causing it to oscillate back and forth. Figure 2a shows the results. Flux tubes are seen 
to oscillate about the equipotential surfaces and in addition they move along it. The 
movement along the equipotential surfaces is due to the tension which pulls the flux 
tube toward the rotation axis. If is sufficiently low, where tension can be totally 
neglected, the flux tube will not move along the equipotential surface, but just oscillate 
both radially and vertically. Figure 2b shows the results for an = 0.1 disk. In addition 
to the tension which acts toward the rotation axis, there is now a Coriolis force which 
acts away from it; these two competing forces determine the fate of the flux tube. For 
example. Mo = 0.3 flux tubes released at = 10 and 14 revert back as soon as they 
cross the dot-dashed curve. This marks the surface where the component of effective 
gravity perpendicular to the rotation axis becomes zero. This is also the surface where 
the component of buoyancy perpendicular to the rotation axis changes sign from being 



14 



away from it outside this surface to being toward it inside this surface. If the field is 
sufficiently large as for the M = 0.7 case, then tension makes the flux tube overshoot 
so much that buoyancy makes it oscillate about the equatorial plane. 



3.2 Radiative and convective disks 



Hereafter, we confine ourselves to the study of fiux tubes which move adiabatically. 
We do not deal with the isothermal case because the generic behaviour of the fiux tubes 
is similar to that in adiabatic or isothermal disks. We choose a radiative thick accretion 
disk (convectively stable) by choosing /3{Wc) = 0.0002 and /3{Wout) = 0.0001 (See § 
2.2). Equation ( ]T7| ) gives the distribution of p in such a disk. Flux rings of Mo = 0.3 
are released from ^ = 5, 7, 10 and 15 and 6 = 89° in a n = disk with zero initial 
velocity. The trajectories are computed by integrating Eqns. (|l])-® assuming that drag 
and accretion are negligible {D = and Uacc = 0) and the fiux tubes move adiabatically 
in the disk (this is ensured by using Eqn. ^). Except for replacing the isothermal disk 
with the radiative disk, the conditions are same as in Fig. 2a. The results are shown 
in Fig. 3a. The fiux tubes oscillate about an equipotential surface as in an isothermal 
disk. For any given entropy gradient there is a maximum Mo which can be confined 
in the disk at a given point. Figure 3a shows that the fiux tube released at ^o = 15 
emerges because the oscillation amplitude increases with ^o (Eqn. ^ given below) and 
the amplitude is so large that it emerges before it can turn back. Here we present a 
simple derivation to show how the oscillation amplitude depends on Mo and ^o- 



Section 3.1.2 described, in detail, how and why the fiux tubes oscillate around 
an equipotential surface, where its density is equal to the external density. (On the 
equipotential surface, densities will be exactly equal if M ~ 0.). The variation of 



density with height is described by Eqn. (|T9D- If a small displacement of the fiux 
tube equalises the external and internal densities, (pt/pe) = 1, we can substitute the 
constants ki and k2 from Eqn. (|19D and approximate Eqn. (|19D in the limit when 

Pe « 1, as 
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Expanding Ke by a Taylor series about the point where the flux tube was released and 
retaining only the first term, we have 



Ke,o + 



OK dW 
dW dx 



0' 



dx ~ 



— w 

1-M 



o 



(25) 



which on simplification gives the amplitude of oscillation as 



dx — 



K^Mo 1 



(26) 



(1 - Mo) K2{dW/dx) ' 



For Mo << 1, the amplitude of oscillations increases with Mo which measures the initial 
departure from hydrostatic equilibrium. Beyond the dot-dashed curve (cf. Fig. 2b), 
{dW/dx) > indicating that the initial displacement of the fiux tube is outward in 
the direction of the pressure gradient force. Similarly, between the inner edge of the 
disk and the center (dW/dx) < 0, hence the initial displacement of the fiux tube is 
toward the rotation axis along the pressure gradient force. The constant K2 = dK/ dW 
indicates whether the disk is adiabatic, convective or radiative. If K2 = 0, the disk is 
adiabatic, and the initial displacement is infinite; in other words, the fiux tube emerges 
out in the direction of the pressure gradient force (Fig. 8 of Paper I). If ^2 > the 
disk is radiative and the displacement, or oscillation amplitude, goes inversely to it. 
The derivation is not valid for K2 < 0. When l3{Wc) is made even slightly greater than 
P{Wout) flux tubes are expelled out as in the convective disks. 

3.3 Inclusion of the effects of shear ampUfication 

We study the behavior of the flux tubes that are anchored in a radiative disk 
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under shear amplification. This is done by amphfying the field through Eqn. ( pOD 
and substituting Eqn. (^I]) for the tension and Eqn. (^) for the buoyancy factor, 
while numerically integrating Eqns. (|l|)-(0). Figure 4 shows the results under similar 
conditions as Fig. 3a; the amplitude of the oscillations increases and also the mean 
surface about which the oscillations take place keeps shifting outward till the fiux tube 
emerges. The process takes several dynamical timescales before the fiux tubes emerge. 

4. IMPLICATIONS OF MAGNETIC ACTIVITY INSIDE THICK 
ACCRETION DISKS 

The magnetic field brought by a fiux tube depends upon the point of release, 
the degree to which the field is saturated in comparison to the equipartition value, as 
well as the location where it surfaces. Because of the uncertainties, we have refrained 
from providing any quantitative description of the field itself. In general, for typical 
values of the parameters, the field is found to range between a few thousand Gauss at 
the inner edge of the disk to a few Gauss at the outer edge. This is quite strong and 
a large number of important observational effects are expected if the chimneys of the 
accretion disks do harbour such fields in the form of loops, coronae and prominences. 
In the case of Blazars (i.e., the optically violent variables and the BL Lac objects), 
the effects would be more prominent since the jets are supposed to be beamed directly 
towards the observers. Indeed, most of the multifrequency observations of these objects 
do invoke the presence of very strong magnetic flares for a satisfactory explanation of 
the magnitude and the time scale of variabilities as well as the energies associated with 
them. Recent observation of TeV gamma rays from a highly variable source MK421 
may also be explained through repeated accelerations in magnetic loops. In the case 
of other radio sources, which are not beamed towards the observer, the presence of 
magnetic flelds in the chimney region can only be inferred through indirect means, by 
comparing the kinematic luminosity of cosmic jets with the mass loss through magnetic 
winds from the chimney surface, for instance. Similarly, if signiflcant magnetic flux 
tubes are present inside the disks, they would play an important role in influencing 
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nuclear compositions both inside the disk, as well as outside. Below, we discuss some 
of these suggestive effects. 

(a) X-ray flickering and optical microvariabilities in AGN: 

Emission from Blazars exhibit high polarization and variability at all wavelengths 
(e.g.. Brown et al. 1989) on time scales varying from roughly months (quiescent emission 
component) at wavelengths longer than 1 cm to weeks at submillimeter wavelengths to 
roughly days (flaring component) in UV and smaller wavelengths. X-rays can also 
vary between days to months. It is believed that the emission originates from a single 
compact region (Gear et al. 1986) of size r ~ 10^^ to lO^^cm and the magnetic field, 
which is largely turbulent (Jones at al. 1985) is S ~ 0.1 — 1 G (Gear et al. 1985, 
Brown et al. 1989). A large number of these objects such as, OJ287, 3C 345, BL 
Lac, and 3C 454.3 etc. show clear evidence of a net magnetic field aligned orthogonal 
to the direction of the jet axis. The spectral index a shows a break at a frequency 
Uf) which is known to evolve with time as the emission mechanism changes. At low 
frequencies, where the loss during the acceleration is not important, a ~ —0.75 but at 
a break frequency, Vf,, when the electrons experience significant radiative loss during 
acceleration, the spectrum steepens to have a ~ —1.25. But if the electrons experience 
a burst of injection and then the radiative loss, the spectrum should have a slope of 
q; ~ — 2 (Khardashev 1962). This variation is indeed observed. For an efficient radiative 
loss which produces steepening at infrared and which is clearly observed, one requires 
the field be very high. All these points to the evidence that a strong magnetic activity, 
very similar to what happens of the solar surface, must be taking place at the chimney 
wall as well. Blandford and Konigl (1979) suggest that if a succession of mild shocks 
continuously reaccelerate the electrons, the flow should behave as roughly isothermal, 
and the field perpendicular to the flow axis at the shock should vary as B oc R~^. This 
condition produces a very high field of around IG at the fiaring region. In a highly 
variable source, such as MK421 (t^,ar ~ 10000s), (Brodie, Bowyer & Tennant 1987) field 
varies as S = Bor'^-^ with B ~ 0.2G. In PKS2155-304, field varies as B ^ B^r'^-^ with 



18 



Bq ~ lOOG close to the core to 70 /xG at Ipc. The trends of these observational results 
agrees with our calculated field values near the chimney. 

In the Introduction of Paper I, we mentioned some models to account for intraday 
variability of the blazars. We propose here that magnetic flares in the chimney can be 
considered to be a major contributor to this variability. Though the non-linear physics 
of the fiares on the solar surface is far from known, it is very reasonable to assume that 
the reconnection of the abundant fiux tubes in the chimney should produce fiares in 
time-scales (tf) of the order of a few times the infall time-scale ti ~ a~^rc/vr, since the 
fiux tubes typically saturate to equipartition value in about ti. Here, a is the viscosity 
parameter, Tc is the center of the disk and Vr ~ (2GM/r)^/^. For a disk around a black 
hole of mass IO^Mq, Tc ~ lOr^ and a ~ 0.01, ~ 1 — 3d. The energy deposited at 
each fiare is about L^Bj/An. Here, L is the length scale of a typical fiare and Bf is the 
average field on the chimney surface. The energy is released in time L/Va, where, Va 
is the Alfven velocity (S|/47rp/)^/^, pf being the density of matter inside the chimney. 
With reasonable estimates of L ~ r^, S ~ lOOG, and p ~ 10~^^gm cm~^, one obtains 
the energy release, time-scale and the rate of release to be given by 2 x lO^^ergs, 0.3d 
and 2 x 10^^ erg s'^, respectively Depending upon the size of the fiares their occurence 
rate will change. Thus a large number of frequent micro- and minifiares are expected 
inside the chimney. Toward the outer edge of the disk, fiares will be very less energetic 
if frequent, or highly energetic but very less frequent. Away from the hole, the fiare 
timescale becomes longer because both the amplification as well as the buoyancy time 
scales should grow as r^/^ and naively the flare time scales should also grow as r^/^. 
(The unknown non-linear physical processes which govern flarings may change this 
dependence significantly.) This property explains the rapid fiickering in X-rays and 
slower variations in optical and radio emissions. 

(b) High energy gamma rays: 

Recently 0.5TeV gamma rays have been detected from a highly active and variable 
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source Mk421 (Punch et al. 1992). It is possible to produce such energetic gamma 
rays, if one assumes that they are originated from first order Fermi accelerated protons 
radiating at the 'thick targets' close to the base of the magnetic field loops, similar to 
what is believed to be the source of hard X-rays and gamma rays on the surface of 
the sun. However, in the chimney region where crowded coronae can be present, the 
particles may be energized in multiple steps, and not just only in two steps as in the 
solar corona (Bai et al. 1983). In this process, moderately energetic protons released in 
a flux tube, traverse back and forth and are repeatedly accelerated at each encounter 
with shock fronts moving parallel to the field lines from each footpoint. 

In the first order Fermi acceleration, the energy gain is proportional to the total 
energy of the particle. Thus protons are more efficiently accelerated than the electrons 
for a given velocity. The energy gain per encounter can be written as (Jokipii 1966): 

/:^W^2(5sllW{(5^^Xp + (5s) 

where j3s is the shock velocity component parallel to the field lines and j3p is the speed 
of the particle in units of velocity of fight, W is the total energy of the particle, jip is the 
cosine of the pitch angle and 7^ is the Lorentz factor of the shock motion. The average 
time interval between two successive collisions with the approaching shock fronts is 
given by (Bai et al. 1983), 

d 

At^— 

cPp < llp> 

where d is the half length of the fiux tube. Assuming shocks move non-relativistically, 
i.e., Ps « /5p! and using the relativistic relation m = mp/{l — PpY^"^, where, nip is the 
rest mass of the proton, one obtains the time variation of total energy during i-th step 
as, 

W^it) - = (WUti? - Wl)exp[A < /.J > {vjd){t - U)] 

where, Wp = rjipC^. Here, Wi-i{ti) is the injection energy of the proton at the begining of 
the i-th step at time ti (i.e., the final energy of the proton at the end of (i — l)-th step). It 
is easy to see that in the non- relativistic limit, assuming the duration of acceleration to 
be equal to the the shock traversal time of the tube d/vg and < /ip >— 0.75, the energy 
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goes up by a factor of exp{3) ~ 20 as in the solar flux tubes (Bai et al. 1983). However, 
as the total energy goes up, this efficiency factor decreases to about exp{1.5) ~ 4.5 
when the ultra-relativistic equation (given above) is used. Assuming initial injection 

energy Eini to be about IMeV (which is just above the energy where the gain rate is 
higher than the Coulomb energy loss rate), a final energy Wn of 1 — 2TeV is obtained 
in n-steps where, 

~ 2WpEiniexp(3n) 

which yields, n ~ 7. Because a large number of steps are required, it may seem to be an 
inefficient process. But the chimney is more confined than an almost flat solar surface, 
and such multiple steps could be possible. Thus, in general, production of high energy 
cosmic rays in the chimney region cannot be ruled out. 

One way to verify if this is indeed what takes place in Mk421, TeV gamma rays can 
be observed after a delay of several days or even weeks after the detection of 1 — lOMeV 
gamma rays, as a considerable time is spent during the acceleration process. (Indeed, 
at each step, d/vs ^ 2GM/vs(? ~ 3d is required, where, d Vg and Vg ~ lO^cm s"^ 
are used.) Furthermore, intermediate energy photons should also be observed. Note 
that the acceleration within shock traversal time does not put any constraint on the 
size of the magnetic flux tube. But a delay between different bands, if observed, would 
constrain the sizes involved. 

(c) Correlations between radio and optical variabilities 

In the case of BL Lac 0716+714, it has been observed that variabilities in radio 
and optical are well correlated (Wagner et al. 1990; Quirrenbach et al. 1991) and some 
kind of quasi periodicities are observed in time scales of 1.2d in the first week and 
of ~ 7d in later weeks of observations. This behavior could be understood using our 
observation of the fact that strength of magnetic fields on the disk surface is well 
correlated with the pressure inside the disk where the flux tube is originated. If, for 
example, the optical variabilities are caused by the formation and disappearance of 
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the hot-spots (non-axisymmetric density perturbations) inside the disk (Chakrabarti & 
Wiita 1993), then the radio flares caused by magnetic fields anchored with these density 
waves should also follow a similar variability. A possible delay, corresponding to the 

buoyancy timescale of a <day (inner disk) to several days (outer disk) is expected. 

(d) Matter supply to cosmic jets via magnetic winds 

Magnetic winds carry almost 10~^^Mq yr~^ from the surface of the Sun. The 
active surface of the chimney of a thick disk, which may have a typical area (~ 2Trrch, h 
being the maxium height of the disk) of ~ lOOOr^ will release roughly ~ AfPfVs 
amount of matter per second, where Ve is the escape velocity of matter. At an average 
distance of r ~ lOr^, and with pf — 10~^^gm cm~^, the rate will be ~ Mq yr~^. This is 
comparable to the typical kinematic fiux observed in a jet. Thus, the magnetic winds 
from the chimney can be responsible for supplying necessary matter in the jet. 

(e) Formation of neutron disks and neutron rich isotopes 

In the presence of a large number of fiux tubes, which are originated deep inside, 
the nuclear composition of the disk could change in buoyancy time scale. In the case of 
ion pressure supported thick disks (Recs et al 1982), the central temperature of the disk 
could be high enough (> 10^°K) that the accreting matter would be photo-dissociated 
completely (Chakrabarti, Jin and Arnett 1987; Hogan and Applegate 1987). However, 
the buoyant fiux tubes will carry with them only the protons and electrons. Neutrons 
left behind can produce a neutron rich disk close to the center. Newly accreting matter 
can capture these neutrons and produce neutron-rich nuclear isotopes which eventually 
escape from the disk by the rising and fiaring fiux tubes. 

5. CONCLUSIONS 

In the present paper, we have been able to show that it is possible to construct 
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physical models of thick accretion disks around black holes which can act as a store 
house of magnetic fields. In isothermal disks or in radiative disks, weaker flux tubes are 
not expelled away. Instead, they continue to oscillate around equipotential surfaces, 
dictated by the forces of gravity and the centrifugal motion, till they are amplified 
by shear and become more buoyant. It is well known that in the case of the sun, 
the flux tubes are stored in the regions near the boundary between the radiative and 
the convective zones, and appear on the surface due to instabilities. Similarly, the 
anchored flux tubes in the disk are also expected to appear close to the chimney due to 
perturbative effects and cause magnetic activities. We have discussed briefly that a large 
number of astrophysical processes, such as supply of matter through magnetic winds, 
variability of Blazars through flaring events, production of highly energetic particles in 
the coronae through Fermi acceleration processes, etc. could be related to the interesting 
magnetic flux tube behavior that we discover. 

In our work, in both Paper I and Paper II, we ignored the effects of reconnection 
of flux tubes within the disk. Reconnection takes place in diffusion timescale, and is 
slower compared to the shear ampliflcation timescale. Therefore, it is not expected 
that it would have any important effect on the dynamics of flux tubes inside the disk. 
However, outside the disk, in the magnetosphere, the Alfvcn velocity is very high (since 
the matter density is low) and thus reconnection is important. Indeed, the flares on the 
solar surface are believed to be due to this reconnection process. 

Though the toy model of the disk we had used is sufficiently reahstic under the 
circumstances mentioned in the text, it is desirable to construct a disk model which 
is fully self-consistent, namely, in which the structure of the disk is determined from a 
minimum number of parameters. For example, we have ignored the effects of nuclear 
reaction and associated energy generation inside a hot disk. In Paper I, we have already 
shown how various parameters that wc use arc closely linked with the 'unknown' viscous 
mechanism operating inside a disk. There is also a large uncertainty in the size and 
fleld distribution of the flux tubes of the accreted matter. Though, once they are inside 
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the disk they are amphfied to the local equipartition value in a dynamical timescale. 
It would therefore appear that a more complete knowledge of the behavior of the flux 
tubes and their influence on the observable phenomena must await an improvement of 

our understanding of the uncertainties mentioned above. On the other hand, since we 
have chosen very reasonable sets of parameters, we believe that the major conclusions 
drawn in our papers (I & II) will remain intact. 



24 



REFERENCES 



Bai, T. et al. 1983, ApJ, 267, 433 
Balbus, S.A. & Hawley, J.F., 1991, ApJ, 376, 214 
Balbus, S.A. & Hawley, J.F., 1992, ApJ, 400, 610 
Balbus, S.A. & Hawley, J.F., 1992, ApJ, 400, 620 
Bisnovatyi-Kogan, G.S. and Blinnikov, 1977, A & A, 59, 111 
Blandford, R.D. & Konigl, A. 1979, ApJ, 232, 34 

Bregman, J.N., 1988, in Supermassive Black Holes, ed. M. Kafatos (Cambridge Uni- 
versity Press), p. 43 

Brodie J., Bowyer S. and Tennant A.: 1987, ApJ, 318, 175 
Brown et al. 1989, ApJ, 340, 129 

Brown, L.M., Robson, E.I., Gear, W.K. & Smith, M.G. 1989, ApJ, 340, 150 

Chakrabarti, S.K. & D'Silva, S. 1994, ApJ, (to appear) 

Chakrabarti, S.K., Jin, L. & Arnett, W.D. 1987, ApJ, 336, 572 

Chakrabarti, S.K. & Wiita, RJ. 1992, ApJ (submitted) 

Choudhuri, A. R. 1990, ApJ, 355, 733 

D'Silva, S. & Choudhuri, A.R., 1993, A&A, 272, 621 

Galeev, A. A., Rosncr, R. & Vaiana, G. S., 1979, ApJ, 229, 318 

Gear, W.K. et al. 1985, ApJ, 291, 511 

Gear, W.K. et al. 1986, ApJ, 304, 295 

Hawley, J. F. & Balbus, S. A. 1991, 376, 223 

Hogan, C. & Applegate, J. 1987, Nat, 330, 236 

Jokipii, J.R., 1966, ApJ, 143, 961 

Jones, T.W. et al. 1985, ApJ, 290, 627 

Khardashev, N.S. 1962, Sov. Astr.- AJ, 6, 317 

Kaisig, M., Tajima, T. & Lovelace, R.V.E., 1992, ApJ, 386, 83 

Lawrence, A., Watson, M.G., Pounds, K.A., & Elvis, M. 1987, Nat, 325, 694 

Landau, L.D. and Lifshitz, F.D. 1959, Fluid Mechanics, Pergamon Press (New York) 

Lightman, A.P. & Eardley, D.M. 1974, ApJ, 187, LI 



25 



Matsumoto, R., Horiuchi, T., Hanawa, T. & Shibata, K., 1990, ApJ, 356, 259 

Moreno-Insertis, F., 1986, A&A, 166, 291 

Paczynski B. 1978, Acta. Astron. 28, 91 

Paczynski, B. & Abramowicz, M. A. 1982, ApJ, 253, 897 

Paczynski B. & Wiita, P. 1980, Astron. Ap, 88, 23 

Parker, E. N., 1955, ApJ, 121, 491 

Parker, E. N., 1979, Cosmical Magnetic Fields, Oxford University Press 
Punch, M. et al. 1992, Nat, 358, 477. 
Quirrenbach, A. et al. 1991, ApJ, 372, L71 

Rees, M.J., Begelman, M.C., Blandford, R.D, & Phinney, E.S, 1982, Nat, 295, 17 

Shakura, N. I. & Sunyaev, R. A., 1973, A&A, 24, 337 

Shibata, K, Tajima, T. & Matsumoto, R., 1990, ApJ, 350, 295 

Wagner, S., Sanchez-Pons, F., Quirrenbach, A., & Witzel, A. 1990, A&A, 235, LI 

Wiita, P.J. 1982, ApJ, 256, 666 



26 



Figure Captions 



Fig. 1 — Trajectories of flux tubes in the i? = ^sin^-2; = ^ cos 9 plane, released at 
{ = 4, 5, 6, 7, 10 and 9 = 89°, with zero initial velocity and no drag, under isothermal 
conditions. The dotted line shows the outer boundary of the disk. The flux tubes have 
Mo = 0.1, the disk has n = 0, and the markers are at time steps of lOOr^/c. (a) The disk 
temperature T^j = 5 x 10^ °K. (b) At Ta = 10^° °K where the tension is the dominent 
force. 

Fig. 2 — Trajectories as in Fig. 1, except that the flux tubes move adiabatically and 
they have Mo = 0.3 (to highhght the oscillations), (a) Here, n — Q i.e., flux tubes do 
not feel any Coriolis force, (b) n — Q.l and flux tubes feel Coriolis force. 

Fig. 3 — Trajectories as in Fig. 2, except that the flux tubes (Mo = 0.1) move adia- 
batically inside a radiative disk; (a) n = 0, (b) n = 0.1. 

Fig. 4 — Trajectories of flux tubes released at ^ = 5 and 7, when the held is shear 
amplifled exponentially. Other parameters are same as in Fig. 3. 
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